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Abstract 

In this paper we consider storage and inventory systems. Our 
aim is to apply and review main results of the fluctuation theory of 
stochastic processes in the context of storage and inventory modeling. 
We describe systems where the inflow is due to a Levy process and 
the outflow is linear and conversely systems where the inflow is lin- 
ear and the outflow is due to a Levy process. For such systems we 
investigate the process of a storage (inventory) level. We give formu- 
las for the probability that the storage level exceeds a certain value. 
This probability is crucial in the management of storage and inventory 
systems. 
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1 Introduction 



Levy processes appear in many theoretical and practical fields where they 
serve as a basic skeleton for a description of certain phenomena. They are 
applied in physics, economics, finance, insurance, queueing systems and other 
branches. Their features like independence and stationarity of increments or 
self-similarity in certain cases for instance permit to apply them to model 
returns of stock prices, claims to insurance companies or an inflow (outflow) 
to the buffer in queueing (telecommunications) systems. Moreover Levy pro- 
cesses serve as a starting point for more complicated models e.g. based on 
stochastic differential equations. The topic of this research paper is Levy 
processes in the context of inventory management. 
Let us recall the notion of a Levy process. 

Definition 1 A stochastic process {X(t),0 <t< 00} with values in TR d is 
a Levy process if 

• X(0) = a.s., 

• X has independent increments, 

• X(t) — X(s) has the same distribution as X(t — s) , that is X has 
stationary increments, 

• X is stochastically continuous. 

The structure of Levy processes is not much complicated. Levy-Ito repre- 
sentation shows their stochastic construction which is the following (see e.g. 
Sato 0) 

X{t) = B{t) + [ x (N t (dx) - tQ(dx)) + [ x N t (dx) + at , 
J\x\<i "'M^i 

where B[t) is a Wiener process, iV is a point process generated by the jumps 
of X that is iV = J2{t:AX(t)^o} ^(t,Ax(t)) • N is a random Poisson measure on 
[0, 00) x {]R d \ 0} with the mean ds x Q(dx) , where Q(dx) is the so-called 
Levy measure on 1R d \ and a e lR d . The abbreviation N t {dx) means 
N([0,t],dx) . 
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In this article a stochastic process describing a storage level in storage and 
inventory systems is investigated. The presented research will distinguish 
between two cases. The first one will assume that the inflow of goods or 
products (we do not distinguish products that is we assume that we store 
one kind of the products) to a storehouse is random due to a Levy process 
and the outflow is deterministic and linear. Such a situation can happen 
when the products peel from a conveyor belt, natural resources are mined 
and then stored or in e.g. fishing industries. A similar situation we can have 
in terminals, dams or telecommunication systems where the inflow is random 
and the outflow is deterministic. This case is particularly interesting as we 
know from the review by Silver 16j that it is uncommon in the scientific 
literature to model supplies (inflow) for storage systems in this manner. 

In the second case we assume conversely that is the inflow is deterministic 
and the outflow is random. Such a situation happens when the customers buy 
a certain product and a storehouse can guarantee deterministic linear supplies 
and the demands of customers are random. In this case we consider the so- 
called inventory systems. From Silver [16( we know that these assumptions 
are quite common in research. 

In this paper we want to review theoretical results which determines the 
probability that the storage level exceeds a certain value. Such a character- 
istic is important in the management of storage and inventory systems. In 
many cases this value can be determined from computer simulation but to 
check the accuracy of the simulation we need theoretical formulas. Moreover 
the theoretical results are interesting from the statistical point of view, that 
is they are necessary to fit a model to real data. They are also noteworthy 
from the mathematical point of view as they are based on deep results from 
the theory of stochastic processes. 

Let us start to describe our models more precisely. In the first case we 
assume that the supply to a storehouse is random and due to a Levy process 
X that is the whole supply to a storehouse until time t is X(t) . From 
the practical point of view we should assume that the processes X is non- 
negative with non-decreasing sample paths but we consider several examples 
where these assumptions will be dropped. Then we can regard that some 
uncertainty (a noise) is added to the outflow from the storehouse. First for 
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the simplicity let us assume that the net supply is the following 

Y(t)=X(t)-t, 

where X is a Levy process with values in [0, oo) and with non- decreasing 
sample paths. Thus the process Z(t) describing the stock level of products 
in the storehouse at time t satisfies the following equation 

Z(t) = Z{0) + Y(t) + f l{Z(s) = 0}ds, (1) 

Jo 

where Z(0) is the stock level at time t — (it can be a random or deter- 
ministic quantity) and 



l{Z(s) = 0} 



1 if Z(s) = , 
if Z(s) > 



which means that the stock level can not be less than zero and if Z(s) = the 
supplies are not greater than demands which are linear. The next proposition 
gives a solution for eq. (JTJ (see e.g. Prabhu 

Proposition 1 The integral equation (TJP has a unique solution in the fol- 
lowing form 

Z(t) = Z(0) + Y(t)+I(t), (2) 

where 

I(t) = [mfY(s) + Z(0)]-. 
Moreover if Z(0) =0 then 

Z(t) = sup[Y(t) - Y(s)} = sup Y(s) . (3) 

s<t s<t 

Proof: Subtracting by the sides the equation ([1]) for times t and r— (the 
left limit at r , r < t ) we get 

Z(t) = Z{t-) + Y(t) - Y (r-) + f l{Z(s) = 0} ds 
> Y (t)-Y(r-), 

where the last inequality follows from the fact that the integral in the equa- 
tion is non-negative. First let us assume that the storehouse is empty at 



some moments until time t (precisely on the time interval (0, t] ). Then we 
can define the last moment when the storehouse is empty (until the moment 
t) that is to = maxjr : r < t, Z(r—) = 0} . Thus to is the last moment 
when the storehouse was empty on the time interval (0, t] and 

Z(t) = Y(t)-Y(t -) 

using the last equation and the fact that / t * l{Z(s) = 0} ds = . Hence by 
the last inequality we obtain 

Z{t) = sup[Y(0 - Y(t-)} (4) 

T<t 

(we put Y(0— ) = 0). Now let us assume that the storehouse is never empty 
on the time interval (0, t] . Then by eq. (JT]) it follows that Z(0) + Y(t— ) > 
for all < r < t hence Z(t) = Z(0) + Y(t) > Y(t) - F(r-) . Thus by © 
we get 

Z(t) = max{sup[y(t) -Y(r-)],Z(0) +Y(t)} 

T<t 

= max{Y{t) -MY(t),Z(0) +Y{t)}, 

T<t 

where in the last equality we use the stochastic continuity of the process Y . 
Thus using the last equality and eq. (pQ) we get 

[ l l{Z(s) = 0}ds = Z(t)-Z(0)-Y(t) 
Jo 

= max{-inf Y(r) - Z(0),0} 

= [MY(t) + Z(0)]-. 

Now let us assume Z(0) = and notice that Yi(r) = Y(j) for r < t (in 
sense of finite dimensional distributions), where Yi(t) = Y(t) — Y(t — r) . 
Then for a given t we have 

Z(t) = Y{t)-MY{s) 

s<t 

= sup[Y (t) - Y (s)} 

s<t 

= sup Yx(t — s) 

s<t 
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= supYi(s) 

s<t 

= supF(s) 

s<t 

where the last equality is in the sense of one dimensional distributions. The 
proof is completed. 

Remark 1 Except the eq. (TJ|) the statements of the above theorem are valid 
for any process X which starts from zero and has non-decreasing trajectories 
and is stochastically continuous that is we can drop the assumption that X 
is a Levy process. 

Remark 2 The equality in |3]) is not in the sense of finite dimensional dis- 
tributions but in the sense of one dimensional distributions. 

Remark 3 The process |IP can be used as a storage level process even though 
the process X has values in 1R and its sample paths are not non- decreasing. 
Then we add some uncertainty to the outflow from the storehouse and the 
process (0j is a solution of a slightly different equation than Thus the 
process |I)) is more general and can serve as a starting point to describe 
storage levels. 

As we mentioned we will investigate a model where the supplies are de- 
scribed by a linear function and the customers demands (outflow) are ran- 
dom. More precisely we assume that the net supplies are the following 

Y{t) = t - X{t) , 

where X is a Levy process with values in [0, oo) and then its sample paths 
are non-decreasing. Hence the process Z(t) describing the stock level at a 
moment t satisfies the following equation 

Z(t) = Z(0) + Y(t) + f l{Z(s) = 0} dX(s) , (5) 

J 

where Z(0) is the stock level at time t — (it can be a deterministic or 
random quantity). 
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Proposition 2 The integral equation |3P has a unique solution of the form 

Z(t) = Z(0) + Y(t) + I{t) , (6) 

where 

I(t) 

Moreover if Z(0) — then 

Z(t) = sup[F (t) - Y(s)} = sup Y(s) . (7) 

s<t s<t 

Proof: The proof goes a similar way as the proof of Prop. HJ 

Remark 4 Except eq. |7p the statements of the above theorem are valid for 
any process X which starts from zero and has non-decreasing trajectories 
and is stochastically continuous, that is we can drop the assumption that X 
is a Levy process. 

Remark 5 The equality in |?P is not in the sense of finite dimensional dis- 
tributions but in the sense of one dimensional distributions. 

Remark 6 The process (0|) can be used as a storage level process even though 
the process X has values in IR and its sample paths are not non- decreasing. 
Then we add a some uncertainty to the inflow to the storehouse and the 
process is a solution of a slightly different equation than Thus the 
process is more general and can serve as a starting point to describe 
storage levels. 

In this paper we are interested in two quantities that is the probability 
of an overflow at time t which is JP(Z(t) > u) and the probability of an 
overflow up to time t which is P(sup s<t Z(s) > u) under assumption that 
the initial storage level is deterministic Z(0) = z > . In several cases we 
present explicit formulas for those probabilities or we give Laplace trans- 
forms of related quantities. In the article D^bicki and Mandjes [5[ a similar 
treatment is conducted but they consider the workload process (in our article 
it is the process Z ) in the stationarity regime. More precisely the random 
variable sup s<00 Y(s) is investigated. They review results which present 



[inf Y( 

s<t 



s) + Z 
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Laplace transforms of the distribution of sup s<00 F(s) for a spectrally pos- 
itive or negative Levy process Y . Moreover the results on the distribution 
of Z stopped at an exponential time are presented and the correlation of 
the workload process Z is derived. In the paper of Debicki and Mandjes 
(lil they also investigate inf s < t Z(s) but as before under assumption that the 
workload process is stationary. Almost all their results are in the language 
of Laplace transforms or double Laplace transforms. Thus our work is dis- 
joint with the work of Debicki and Mandjes Q| and we emphasis on explicit 
formulas. 



2 Probability of an overflow in storage sys- 
tems 

In this section we will investigate the process Z from Prop. [TJ Thus un- 
der assumption Z(0) = (see eq. ([3])) we are interested in the following 
probability 

W{Z{t) >u) = P(sup(X(s) - s) > u) 

s<t 

which is the probability that the stock level exceeds u at time t . First we 
assume that the process X is spectrally positive, i.e. its Levy measure Q 
has the support on the positive half line (this means that its jumps are only 
positive) and has a finite variation which is equivalent for the Levy process X 
that Jq 1 xQ(dx) < 00 . The next theorem is a direct consequence of Michna 



9| and lCj which is originally due to Takacs 19j in a little more general 
setting (for the processes with non-negative interchangeable increments that 
is in our case X has to have the finite variation). 

Theorem 1 Let u > , Z(0) — and X be a spectrally positive Levy 
process with the finite variation. Then 

P{Z(t) > u) = P(sup(X(s) - s) > u) 

s<t 

= ¥{X{t)-t> u) + 

H f(u + S,S^ rt-s 



Jo t — s Jo 
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where f(x,s) is the density function of X(s) provided that it exists and 

P(Z(t) > 0) = 1 - - / P(X (t) <x)dx. 
t Jo 

Example 1 Let X be gamma process. Then its Levy measure is the follow- 
ing 

Q(dx) = ( x > 0} f—?L \ olx , 
x \ b J 

where a > and b > . This process is a subordinator e.g. its trajectories 
are non- decreasing. The density function of X(s) is of the form 



f( x , s) = - — . x as ~ l exp -- l{x > 0} . 
v ' b as T(as) V bj 

Thus the quantity of the inflow until time s is distributed due to the above 
density function. Hence using Th. U\we get for u > (see also Dickson and 
Waters %) 

V{Z{t) > u) 
= ¥(X(t) -t> u) 

+ /' P(X (t -s) <t- s)f(u + s, s) ds 
Jo 

-ah ( P(X(t -s + l/a) <t - s)f(u + s, s) ds , 
Jo 



where 



Example 2 Let X be inverse Gaussian Levy process. Then its Levy mea- 
sure is the following 

Q(dx) = exp(- j-73;) dx 

V 27TX 6 * 

where 5 > and 7 > . This process is a subordinator. The density 
function of X(s) is the following 

s ) = ^=L" 3 / 2 exp f -I^^x" 1 + 7 2 x)l l{x > 0} . 
V27T 12 J 



As in the previous example for u > we obtain 



P(Z(t) > u 

5te^ 5t 



[ x~ 3/2 exp (--(^VaT 1 +7 2 x)l dx 
Ju+t 12 J 



y/2n Ju+t 
S 2 e^ 5t ^ , f 1 



J s(u + sy 3/2 exp j--(5V(it + s)~ l +7 2 (m + s))| ds 



2tt 

• £~\(t - s)x- 3/2 - x~ l/2 ) exp {-\(S 2 (t - s) 2 x- 1 + 7 2 x)| dx . 
Similarly one can consider generalized inverse Gaussian Levy process. 

Now let us investigate the case when X is a spectrally positive Levy 
process with the infinite variation which means that J xQ(dx) = oo . Let 
us notice that in this case the trajectories of X are not non- decreasing which 
means that in our model we add a some uncertainty to the outflow. The next 
result is due to Michna 191 and 10 . 



Theorem 2 Let u > , Z(0) = and X be a spectrally positive Levy 
process with the infinite variation. Then 

P(Z(t) > u) = P(sup(X(s) - s) > u) 

8<t 

= P(X(t) -t> u) + 

r ds r s &(x(t -s)<x)dx, 

Jo t — S J-oo 

where f(x,s) is the density function of X(s) provided that it exists and 

X*(Z(t) > 0) = 1. 
For the proof see Michna 0] and (Io| . 

Example 3 Let X be an a -stable spectrally positive Levy process with 1 < 
a < 2 . Its Levy measure is of the form 

Q( dx ) = ° } dx 

where a > . The density function of X(s) exists but it can not be presented 
as an elementary function (see e.g. Zolotareav 12W or Nolan I IB])- There 
exist integral or series forms of f(x, s) so applying Th. [H one can compute 
numerically probability that the stock level will be greater than u at time t . 
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Example 4 Let X be a tempered a -stable spectrally positive Levy process 
with 1 < a < 2 that is its Levy measure is the following (see Koponen f^J) 

where a > and A > . T/ie density function of X(s) one can determine 
by the inverse Fourier transform of the characteristic function. Thus by Th. 
IE we can compute numerically probability that the stock level will be greater 
than u at time t . 

According to Prop. [I] the stock level at time t for Z(0) = z (z > is 
deterministic) is the following 

Z(t) = z + X(t) - t + [inf (X(s) -s) + z\- . 

s<t 

Let us put 

Ee -™m) = e ^(«) ; ( 8 ) 

where w > and 

t{u) = mf{t > : Z(t) > u] . 

Let us notice that t{u) is the first epoch when the stock level exceeds u . 
Moreover probability that the stock level exceeds u until time t is 

P(sup Z(s) >u) = P(t(m) < t) . 

8<t 



Next theorem is due to Pistorius (see Pistorius [13| see also Avram et al. 



Theorem 3 For < z < u 

Ee -rr(u) = K {r)( u _ z }_ rW (r)( u _ z )W {r \u)/W {r) ' + {u) 

where r > and 

K {r \x) = 1 + r f X W {r \y)dy, 



o 

where W^ r \y) is the so-called r -scale function on [0, oo) determined by 
Laplace transform 

°° e-"WW(y) dy = 

o ip{w) — r 

11 



defined for w > $(r) , where $(•) is the right inverse function to ip(-) (see 
eq. (EPJ and W { + ] {u) is the right derivative of W^ r '{x) at x = u . 

The above theorem gives the full solution to the problem of the overflow 
of the stock level in the model with the deterministic linear outflow. The 
disadvantage of the above formula is that we can obtain the probability only 
numerically by inverting Laplace transform of complicated functions. 

From the above theorem one can easily find the expected value of the first 
epoch when the store level exceeds u . Let us put W(x) = W^\x) . 



Corollary 2.1 



W(u) 

JEr(u) = W(u-z) w ^--W(u 



where W + {u) is the right derivative of W(x) at x = u and W(x) = 

3 Probability of an overflow in inventory sys- 
tems 

In this section we will assume that the supplies are linear and the outflow 
is a spectrally positive Levy process (its jumps are only positive) that is the 
net supplies until a moment t are the following 

Y(t) = t - X(t) 

where X is a spectrally positive Levy process. Let us notice that Y is now 
a spectrally negative Levy process. From the mathematical and practical 
point of view the situation is different than in the model with the linear 
outflow. Under these assumptions we give a complete solution for any initial 
stock level. This is highly desirable as it gives us a general model, for any 
outflow modeled as a Levy process. Furthermore this is also interesting 
as this probability is directly linked to the term of a service level within 
inventory management see Silver et al. [13] ■ The service level corresponds to 
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the percentage of customer demand (in the form of volume) serviced directly 
from stock. In the simplified case when Z(0) = the following probability 



W{Z{t) >u) = P(sup(s - X(s)) > u) 



can be derived from Kendall's identity which is valid for spectrally negative 
Levy processes, here Y(t) = t — X(t) is spectrally negative (see Kendall 0] 
or e.g. Borovkov and Burq and the references therein, see also Takacs 



19l | for the case when X has the finite variation). For an elementary proof 
of an explicit formula for the above probability in the case of the spectrally 
negative a -stable Levy process without any linear drift with 1 < a < 2 see 



Michna In a more general setting we find the probability that the stock 



level will exceed u until the moment t that is 

P(supZ(s) > u) 



s<t 



if Z(0) = z for < z < u and z is deterministic. However as before the 
solution will not be in an explicit form but as Laplace transformations of 
scale functions. 

According to Prop. [2] the stock level at time t is the following 
Z(t) = z + t- X(t) + [inf(s - X(s)) + z\- . 

s<t 

Let us put 

Ee »ni) = e ^M ? ( 9 ) 

where w > and 

t{u) = mf{t > : Z(t) > u] . 

Let us notice that t(u) is the first epoch when the stock level exceeds u . 
Moreover probability that the stock level exceeds u until time t is 

P(sup Z(s) >u) = P(t(m) < t) . 

s<t 



Next theorem is due to Pistorius (see Pistorius |13| ) 
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Theorem 4 For < z < u 



K<r)(u) 



K^ r \x) = 1 + r / W (r) (y)cfy 



where r > and 



where W^'\y) is the so-called r -scale function on [0, oo) determined by 
Laplace transform 

°° e - wy W (r \y)dy ' 



^(w) — r 

defined for w > $(r) 7 where $(•) zs i/ie n'a/ai inverse function to ip(-) (see 
eg. W). 

Remark 7 eg. (^) and (0|) i/ie function if) and the scale functions in 
Th. and |7] are the same. 

Similarly as in Th. [3] the above formula permits to determine the probability 
of the overflow of the storehouse but we do not avoid numerical calculations. 
For some Levy processes the scale function is available in an explicit form 
that is for compound Poisson processes, a -stable processes with a e (0,2) , 
tempered a -stable processes and Brownian motion (see e.g. |2j and the 
references therein). Generally one can not obtain an explicit form of scale 
functions. However there are numerical methods which permit to determine 
the scale function of a general spectrally negative Levy process (see Surya 



18( | for Laplace inverse algorithm). 

From the above theorem one can easily find the expected value of the first 
epoch when the store level exceeds u . 

Corollary 3.1 

Er(ii) = W(u) - W(z) . 
where W(x) = W^(x) and W(x) = J* W(y) dy . 

An interesting aspect of this research is the relative ease with which asym- 
metrical distributions of the outflow can be incorporated. This is interesting 



since we know from e.g. Bobko and Whybark [3| and Zotteri [2JJ that the 
asymmetrical distributions are particularly difficult to manage and tend to 
lead to high costs of the systems. 
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